Microscopically derived Ginzburg-Landau theory for magnetic order in the iron 
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We examine the competition of the observed stripe spin density wave (SDW) with other commen- 
surate and incommensurate SDW phases in a two-band model of the pnictides. Starting from this 
microscopic model, we rigorously derive an expansion of the free energy in terms of the different 
order parameters at the mean-field level. We show that three distinct commensurate SDW states 
are possible and study their appearance as a function of the doping and the electronic structure. 
We show that the stripe phase is generally present, but its extent in the phase diagram depends 
strongly upon the number of hole Fermi pockets that are nested with the electron Fermi pockets. 
Electron pockets competing for the same portion of a hole pocket play a crucial role. We discuss 
the relevance of our results for the antiferromagnetism of the pnictides. 

PACS numbers: 75.30.Fv, 75.10.Lp, 74.70.Xa 
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I. INTRODUCTION 

The superconductivity of the iron pnictides contin- 
ues to fascinate the condensed matter community^— Be- 
cause of their high critical temperatures, particular inter- 
est has focused upon the so-called 1111 family 4 - i?FeAsO 
and the 122 family 5 - AFe2As2 (R and A are rare-earth 
and alkaline-earth elements, respectively), which become 
superconducting by chemical doping or under pressure. 
The parent compounds are antiferromagnetsp 3 - with 
stripe-like magnetic order with respect to the lattice of 
Fe atoms. Furthermore, the antiferromagnetism is inti- 
mately linked to an orthorhombic distortion of the crys- 
tal,—^ as evidenced by the coincidence of the ferromag- 
netic direction with the shorter crystallographic axis in 
all 1111 and 122 parent compounds. It has been argued 
that the same condition that favors stable stripe order 
also implies a nematic transition above the Neel temper- 
ature Tjy, where the magnetic fluctuations on each sub- 
lattice become locked into a stripe configuration^r— and 
which produces the orthorhombic distortion via the mag- 
netoelastic coupling^ - - — The mechanism for stabilizing 
the stripe order is therefore a key problem in pnictidc 
physics. 

The microscopic origin of the antiferromagnetism in 
the pnictides has been approached in a number of dif- 
ferent ways. Frustrated local moment models for the Fe 
spins can reproduce the observed magnetic order) 2 ' 13 : 14 
but the evidence for the metalicityi^— and relatively 
weak correlations 1 ^ of the parent compounds, and the de- 
velopment of incommensurate (IC) magnetic order upon 
doping^ suggest an itinerant description. Ah initio cal- 
culations predict ) 20 : 21 and angle-resolved photocmission 
and magneto-oscillation experiments confirm ; 22 : 23 that 
the Fermi surface of the pnictide parent compounds have 
quasi-two-dimensional nested electron and hole pockets. 
Such a system is known to undergo an excitonic insta- 
bility towards a spin-density-wave (SDW) state , 24 ' 25 as 
for example in chromium i 26 ' 27 Many authors have in ad- 
dition emphasized the importance of the complicated or- 



bital structure of the Fermi surfaces ; 28-36 but key aspects 
of the physics are nevertheless well understood on the ba- 
sis of simpler orbitally trivial "excitonic" models i 31 ' 37 ~ — 



Most itinerant models of the pnictides display at least 
two nesting instabilities at different wavevectors. There 
is hence competition between the stripe magnetic order 
and other SDW phases. Within a minimal two-orbital 
model,— it has been shown that doping away from half- 
fillin g 29 ' 30 or a relatively large ratio of the Hund's rule 
coupling to the Coulomb repulsion^ - can stabilize the 
stripe state. For excitonic models, on the other hand, 
Ercmin and Chubuko\^ 3 - have demonstrated that the el- 
lipticity of the electron pockets or interactions between 
the electron bands can stabilize the observed SDW state. 
The stripe order was nevertheless found to be rather sen- 
sitive to the number of Fermi pockets involved in the 
SDW, and its extent in the phase diagram remains un- 
certain. Competition with a different excitonic instability 
has also been proposed to stabilize a stripe SDW 4^ 



In this paper we present a systematic study of the 
magnetic order in the popular two-band excitonic model 
of the pnictidesi 31 ' 32 ' 41 ' 42 ' 45 ' 50 Keeping only interaction 
terms which lead to an excitonic state, we construct the 
Dyson equation for the single-particle Green's functions 
in an arbitrary commensurate SDW phase treated at the 
mean-field level. By iterating the Dyson equation, we 
obtain approximate forms for the self-consistency equa- 
tions for the order parameters valid near Tjv, which we 
then integrate to obtain a Ginzburg-Landau expansion 
of the free energy. From this we determine the phase 
diagram for several choices of the non-interacting band 
structure, and show that three different commensurate 
SDW states are possible. We conclude with a discussion 
of the relevance to the magnetism of the pnictidc parent 
compounds. 
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FIG. 1: (color online) (a) Dispersion of the electron and hole 
bands for the cases of a single hole pocket at the V point 
(fh = 0) and of hole pockets at both the F and M points 
(£h = 0.975). Panels (b) and (c) show the corresponding 
Fermi surfaces and the nesting vectors Qi = (7r/a,0) and 
Q2 = (0, ir/a). In all panels we set fi = 0. 



II. MODEL 

We model the FcAs planes as a two-dimensional inter- 
acting two-band system, where one band has electron-like 
Fermi pockets, while the other has hole-like Fermi pock- 
ets. Including only interaction terms which directly lead 
to an excitonic instability, we write the Hamiltonian as 

H = E{( e k-A0cL^ + (4-M)/L/ kCT } 



k.a 



C k+q,<T C k(r/k'-q,cr'/k'(T' 

cr.er' 

{ C k+q.t C k'-q,i/k',i/k,t + H - C '} .C 1 ) 
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where ci {ft*-) creates a spin-cr electron with momen- 
tum k in the hole-like (electron-like) band. In terms 
of the single-Fc unit cell, we assume the dispersions 
e£ = e c + 2t c (l - £ fe )[cos(& x a) + cos(k y a)} + 2t c £ h [l + 
cos(k x a) cos(k y a)] and e£ = ey + tf t i cos(k x a) cos(k y a) — 
tf,2£,e[cos(k x a) + cos(kya)], where a is the Fe-Fe bond 
length. In units of t c , we take e c = —3.5, £/ = 3.0, 
tf i = 4.0, and t/2 = 1.0. We plot representative band 
structures and Fermi surfaces in Fig. [TJ The dimension- 
less quantities £ e and £/i are key control parameters: £ e 
controls the cllipticity of the electron pockets, while vary- 
ing fa from to 1 tunes the band structure from a sys- 



tem with a single hole pocket at the T point to a system 
with equally large hole pockets at both the T and the M 
points. For £h ~ each electron pocket is nested with the 
hole pocket by only one of the orthogonal wavevectors 
Qi = (n/a,0) and Q 2 = (0,tt/o) [Fig. [[]»], while for 
£h « 1 both electron pockets can nest with a hole pocket 
at each nesting vector [Fig. QJc)] . A system with a single 
hole pocket and two electron pockets was proposed as a 
minimal model of the pnictides in Ref. |43| . and has been 
examined by a number of authors . 42 ] 44 i 47 ] 50 On the other 
hand, a Fermi surface with hole pockets at the T and 
M points is realized in the minimal two- orbital model of 
the pnictides 2 ^ and this situation has been extensively 
studied.— ~ 32 ' 43 ' 47 Furthermore, it is also of relevance to 
more sophisticated orbital models where in addition to 
the d x z / d y z-derived hole pockets at the T point there is 
usually also a d^y-derived hole pocket at the M point , 33 ' 34 
which may play an important role in generating the SDW 
order i 35 ' 36 Although the orbital content of the T and M 
hole pockets are very different, mean- field studies sug- 
gest that the SDW instability is primarily determined by 
the nesting properties ) 35 ' 36 hence justifying the orbitally- 
trivial excitonic model used here. 

Equation (fTJ) contains a density-density interaction 
and a term describing correlated transitions between 
the electron and hole bands, with contact potentials 
gi and g 2 , respectively. At sufficiently low tempera- 
tures, the system is unstable against an excitonic SDW 
with effective coupling 3sdw = 9i + 92 > i 25 ' 37 ' 42 
For our system the excitonic SDW state has two order 
parameters corresponding to electron-hole pairing with 
relative wavevector equal to Qi and Q2, i.e., Ai = 



Ea,^ 

A 2 = 



a,0 



C k+ Ql ,a/k/3) and 



Ea,/3 ' 



^2,a,/3 = (V^)EkEa,/3^«.0( C k+Q 2 ,a/k/3)! 

where a is the vector of the Pauli matrices. Ai and 
A2 are related to the magnetization of each Fc sub- 
lattice by m a = Ai + A 2 , mj, = Ai — A 2 . When 
both Ai and A 2 are non-zero, therefore, the magne- 
tization is the superposition of two SDW states with 
orthogonal ordering vectors. It has been pointed out 
that in the case that Ai • A 2 7^ one has to in- 
troduce additional charge-density-wave (CDW) order 

and Sf = 

Q 2 ffi 



,t 



parameters S c = (1/V) E^^+q,,^/ 

0-/V) E*,Ml+Q,,J*r)> whcre Q3 = Ql + 



III. FREE ENERGY EXPANSION 

We define the single-particle Green's functions of the 
excitonic SDW system by 

G^A^n) = ! ' dr(T T a k+Q Jr)bl a ,(0))e^ T , 
Jo 

(2) 

where a, b — c, /. Treating the SDW and CDW orders at 
the mean-field level, we write the Dyson equation for the 
Green's functions as 
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-9sdw E E A ™^^ Ga(0) ( k + Q»' lw ™) G Q'+Q m ^,-'( k ' lw »)' ( 3 ) 

m— 1,2 a 

i 



where c — f and / = c, and the Green's functions of the 
non-interacting system are G a ^ (k, iu n ) = (iuj n — e£ + 
/z) . The order parameters can be expressed in terms 
of the Green's functions as 

A m , CT)CT , = ^E^E^v( k >^) eiWn0+ > ( 4a ) 

k " n 

8 V=CJ = ^E^E^.r,-( k '^) e ^ 0+ - ( 4b ) 

k.er n 

where /3 = l/fc^T. In general, it is not possible to analyt- 
ically solve Eq. for the Green's functions. By iterating 
the Dyson equation, however, we are able to expand the 
Green's functions in terms of the order parameters. In- 
serting this expansion into the self-consistency ( "gap" ) 
equations and truncating it above a given order, we 
hence obtain an approximate form of the self-consistency 
equations valid close to Tjy (assuming a second-order 
transition, as is the case here). Since the self-consistency 



equations are obtained from the stationary points of the 
free energy with respect to the order parameters, we can 
construct a Ginzburg-Landau expansion for the free en- 
ergy F by integrating them, 

F = F + a(|A 1 | 2 + |A 2 | 2 )+/3 (|A 1 | 4 + |A 2 | 4 ) 
+A|A 1 | 2 |A 2 | 2 + /3 2 |A 1 .A 2 | 2 
+ (7c<5 c + 7/#/)Ai • A 2 + a c fS c Sf + a c 5 2 c + a/<5 2 , 

(5) 

where Fq is independent of the order parameters. Wc 
keep only second-order terms involving the CDW order 
parameters, since the system is far away from a pure 
CDW instability and a CDW only emerges as a secondary 
order parameter^ We also neglect gradient terms since 
we are only interested in homogeneous states. The co- 
efficients in Eq. ([3]) are written in terms of the non- 
interacting Green's functions as follows: 



ft 



7c 
7/ 

a c 



1 + ^ E ^ E GC(0) ( k > ^)° m ( k + Qi. 

k n 



2gsDw 

E \ E [ GC(0) ( k ' ^«) G/(0) (k + Qi,«0 



SDW 



E \ E \ GC(0) ( k ' ^™)G c(0) (k + Q 3 , iu n ) Gf<® (k + Qi , uu n ) 



V 



+ G / ^(k + Qi, iu n )GtW(k + Q 2 ,iu n ) [G c ^(k, iw„)] 

- G c ^(k 1 icu n )G<°\k + Q 3l icu n )G f ^{k + Q 1 ,tuj n )G f(0 \k + Q 2) iw«)} , 

%^ E 7? E GC(0) ( k ' ^«) GC(0) ( k + Q3, ™ n )G m (k + Qi , ( °) (k + Q 2 ^„) 



^fpw £ 1 £ G c ( o )(kj , Wn)G /(o) (k + Ql> ia , n)G /(o) (k + q 2) lWn) ; 



E 4 E GC(0) ( k ' ^«) GC(0) ( k + Q3, iw„)G'<°>(k + Ch, iu n ) , 



5i 



4 S ^ E G/(0) ( k + Qi- -«) G/(0) (k + Q 2 , «*.) 

k n 

GC(0) ( k ' iw -) GC(0) ( k + Q3, ™») ■ 

k n 



(6a) 
(6b) 



(6c) 
(6d) 

(6e) 

(6f) 

(6g) 
(6h) 

(6i) 
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The CDW order parameters can be integrated out, re- 
sulting in the renormalization of 



h = h + 



a 



cf 



4a c a f 



(7) 



The Ginzburg-Landau expansion of the free energy, 
Eq. ([5]), and the expressions for the coefficients in terms 
of a specific microscopic model, Eq. ([6]) , are the first ma- 
jor results of our paper. 



IV. PHASE DIAGRAM 

The free energy in Eq. §5§ admits three possible com- 
mensurate SDW states which we name following Ref. [29|: 

• A magnetic stripe (MS) state where only one of the 
excitonic order parameters is nonzero, e.g., Ai =/= 0, 
A2 = 0. This corresponds to the ordering in the 
pnictides. This state minimizes the free energy if 
2/3 <min{/? 1 +^ 2 ,/3 1 }. 

• An orthomagnetic (OM) state where |Ai| = |A 2 | 
and Ai _L A 2 . This corresponds to a "flux" type 
ordering of the magnetic moments. This state min- 
imizes the free energy if /3i < mm{2/3 , /?i + /3 2 }- 

• A spin and charge order (SCO) state where |Ai| = 

A 2 | and Ai|| ± A 2 . In this state only one sublat- 
tice of the Fe plane has non-zero moments, which 
order in a checkerboard pattern. When gi / 0, 
the spin order induces a CDW with ordering vec- 
tor Q3. This state minimizes the free energy if 

fa +02 <min{2A,,M. 



From close examination of Eqs. (|6b[) and (|Bc|) we observe 

that if ej. = £k + q 3 or e£ = £k+Q 3 * or a ^ k we have 2/?o = 
and hence the MS and OM states are degenerate. 
These conditions are satisfied for the electron and hole 
bands at £ e = and £h = 1, respectively. In particular, 
if £h 7^ 1 the degeneracy of the MS and OM states is 
lifted by arbitrarily small ellipticity of the electron Fermi 
pockets, as pointed out in Ref. l43l 

The free energy in Eq. (|S|) allows us to determine the 
phase diagram of the model close to TV. In Fig. [2] we 
present phase diagrams showing the ordered state re- 
alized at a temperature T = infinitesimally below 
Tjy as a function of £ e and of the doping relative to 
half- filling, dn, for various values of In construct- 
ing the phase diagrams, we adjust (?sdw such that for 
each value of £ e the maximum critical temperature of 
the commensurate SDW as a function of the doping Sn 
is k B T^ pt = 0.0646 t c , which gives a reasonable ratio 
of fcsT^ p to the bandwidth. The variation of the op- 
timal doping level 5n opt , where TV is maximal, with 
£ e is shown by black dotted lines. The boundaries be- 
tween the different commensurate SDW phases are deter- 
mined by the conditions on the /3q, j3\ and /3 2 mentioned 
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FIG. 2: (color online) Magnetic order at T = T N as a function 
of £ e and 8n for (a) £ h = 0, (b) & = 0.95, and (c) £ h = 0.975. 
Magnetic phases are as denned in the text. At £ e = the 
thick solid line indicates degenerate MS and OM solutions. 
For £ e 7^ 0, solid lines indicate phase boundaries, while the 
dotted line indicates the optimal doping Sn opt . 
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S. = 1,S»=0 5, = 1,5,= 0.95 ^=1,^ = 0.975 




(a) 5)1 (b) 5n (c) 8n 



FIG. 3: (color online) Critical temperature Tjv of the SDW states as a function of doping Sn for £ e = 1 and (a) ^ = 0, 
(b) £h = 0.95, and (c) £h = 0.975. The temperature is scaled by the maximum temperature for commensurate SDW order 
k B T° pt = 0.0646 t c . 



above, where the coefficients in Eq. ((BJ) were evaluated 
for gi = g 2 and using a 1000 x 1000 k-point mesh. In 
all phase diagrams we find regions where IC SDW or- 
der occurs. Since the IC SDW ordering vector is likely 
close to the commensurate SDW vector^ the bound- 
ary between the two phases is determined by solving 
1 + (gsow/V) Ek(V/3) E„ G c(0) (k, tw n )G'<°> ( k + Q x + 
<5q, iw n ) = for the critical temperature of the Qi + <5q 
SDW state, where <5q = (O.OItt/o, 0), (0, O.OItt/o). When 
the critical temperature of the Qi +<5q SDW exceeds that 
of the commensurate SDW, an IC SDW is assumed to be 
realized. We similarly find the critical doping for which 
there is no IC SDW order, and the system remains para- 
magnetic (PM) down to zero temperature. Note that we 
disregard states with T/v < 0.05 T^ p . 

We first consider the phase diagram for £/,. = 
[Fig. HJa)], which corresponds to a Fermi surface with 
a single hole pocket and two electron pockets at optimal 
doping as shown in Fig. Hfb). A commensurate SDW 
state is realized here for Sn « Sn opt ±0.025 for all £ e . At 
£ e = the condition e£ = Ck+Q 3 ^ s satisfied, and hence 
the OM and MS states are degenerate. These states have 
the lowest free energy at undcrdoping and near optimal 
doping, but at overdoping the SCO is realized. Upon 
switching on a finite £ e , the degeneracy of the MS and 
OM states is lifted, and the MS state is found to have the 
lower free energy near optimal doping and at overdoping, 
while the OM state is stable at undcrdoping. The SCO 
state is rapidly suppressed by a finite £ e . 

The phase diagrams for the case of two hole and two 
electron Fermi pockets [see Fig. HJc)] is shown in Figs. 
d;b) and (c) for £ h = 0.95 and £ft = 0.975, which cor- 
respond to hole pockets of lesser and greater similar- 
ity, respectively. We note that the interaction strength 
Ssdw needed to produce the SDW state is roughly a 
third smaller than for the single-hole-pocket scenario. 
For small £ e , a commensurate SDW is nevertheless real- 
ized over a much greater doping range than in the single- 
hole-pocket case. The OM phase is stable near optimal 
doping, with the MS phase found at moderate doping, 
and the SCO found at stronger doping. At larger val- 



ues of £ e , however, we find a strong tendency towards IC 
order in the ^ = 0.95 case, with commensurate order 
completely absent for £ e > 1.75. In contrast, the com- 
mensurate SDW in the ^ = 0.975 case is present for all 
£ e and is always realized about optimal doping. 

In Fig. [3] we plot the critical temperature of the SDW 
states as a function of doping Sn for constant-^ cuts 
through the three phase diagrams in Fig. [3J For the 
£h = case [Fig. [3[a)] we note that there are substantial 
IC SDW "shoulders" to the commensurate SDW dome, 
which extend up to T « 0.75 T^ pt . Although IC SDW 
states are also found at strong undcrdoping or overdop- 
ing in the ^ = 0.975 case [Fig. [3^c)], they are realized 
over a smaller doping range relative to the commensu- 
rate states and do not extend to such high temperatures 
compared to the singlc-holc-pocket scenario. As shown 
in Fig. [3jb), however, slightly reducing ^ leads to IC 
states appearing at optimal doping. 



V. DISCUSSION 

To summarize our main results, we have shown that 
in a two-band model of the pnictides there are three dis- 
tinct commensurate SDW states: The MS, OM, and SCO 
phases. In a model with a single hole pocket and two elec- 
tron pockets, the MS state dominates the phase diagram, 
but the OM and SCO phases are possible away from op- 
timal doping. For a model with two hole pockets, in con- 
trast, the OM phase is stable at optimal doping, although 
the MS phase remains at under- and overdoping. Since 
only the MS state is observed experimentally, we hence 
conclude that the model with a single hole pocket gives 
a more reasonable description of the physics. We never- 
theless note that such a model displays a rather strong 
tendency towards IC SDW order, which is not observed 
experimentally. 19 

We consider the results for the single-hole-pocket case 
in more detail. In agreement with Rcf. |43| we found 
that MS order was realized near optimal doping for ar- 
bitrarily small cllipticity of the electron Fermi pockets. 
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underdoping optimal doping overdoping 




FIG. 4: (color online) Schematic diagram of the nesting of the 
two electron pockets (blue dotted and red dashed lines) with 
a single hole pocket (black solid line). We show the situation 
for (a) underdoping, (b) optimal doping, and (c) overdoping. 
The small shaded circles indicate the region of best nesting of 
the electron pockets with the hole pocket. 



Away from optimal doping, however, states consisting of 
a superposition of commensurate SDWs with orthogo- 
nal ordering vectors Qi, Q2 are realized. This can be 
understood via the following argument. At strong un- 
derdoping, we expect that the best nesting between the 
hole pocket and the elliptical electron pockets occurs for 
the states near the major axis of the electron pockets, 
as shown in Fig. HJa) ■ Similarly, for strong overdoping 
the best nesting occurs for the states near the minor axis 
of the electron pockets [Fig. H[c)]. In both cases, the 
SDW gaps for the two nesting vectors involve states far 
apart on the hole Fermi surface, and so there should be 
little competition between them. Near optimal doping, 
however, the nested electron Fermi pockets compete for 
the same states on the hole Fermi surface [Fig.^b)], and 
hence it is more favorable for only a single electron pocket 
to participate in the SDW. We note that the variation of 
the nesting "hotspots" with doping is expected to have 
significant consequences for the a-b resistivity anomaly 
in the pnictidesi^ 3 - 

The addition of a second hole pocket at the M point 
strongly affects the physics. The doping range of com- 
mensurate SDW states is significantly expanded, and the 
OM state is realized near optimal doping with a MS phase 
appearing upon doping. This is consistent with previous 
investigations of the two-orbital modeL 29 ' 30 In contrast, 
our results are inconsistent with the degenerate OM and 
MS phases found in Ref.H^. This is due to the fact that in 
the model of Ref . the T hole pocket is mapped exactly 
onto the M pocket by translation of Q3, i.e., the degener- 
acy condition e£ = £k+Q 3 ^ s always satisfied. Despite the 
apparent differences, the phase diagram can be under- 
stood in a similar way to the single-hole-pockct case. At 
optimal doping the nesting of the electron pockets is not 
optimal for either hole pocket, but instead corresponds 
to overdoping for the smaller pocket and underdoping 
for the larger pocket, thus allowing the OM phase. In- 
deed, this perfectly describes the nesting properties of 
the two-orbital model at half-filling^ Upon doping the 
system, the nesting with one of the hole pockets is opti- 
mized, while the nesting is the other hole pocket becomes 



extremely poor. As such, only a single hole pocket par- 
ticipates in the nesting, and so the MS state is stable. 
The relative size of the two hole pockets is thus crucial: 
If the two pockets are too dissimilar in size, there will be 
a tendency for the T[f(Sn) curves in Fig. [3] to split into 
two separate domes with commensurate order near their 
maxima. The strong tendency towards an IC SDW state 
at optimal doping suggests that the £/, = 0.95 case is close 
to this limit. For somewhat lower i.e., when the hole 
pocket at the M point is much smaller than the pocket at 
the T point, the pocket at the M point only shows good 
nesting with the electron pockets at very large hole dop- 
ing. At realistic doping levels, the smaller hole pocket 
is essentially irrelevant for the SDW formation and the 
singlc-holc-pockct model is applicable. 

Finally, we consider the implications of our results for 
the hypothesized nematic state in the pnictides. Indeed, 
a major motivation for our study is the connection of this 
phase with the MS SDW stated— This can be seen via 
the following naive argument: After integrating out the 
CDW, we write the free energy Eq. ([5]) in terms of the 
sublattice magnetizations m a and mj, 

F = F + ia(|m a | 2 + |m h | 2 ) 

+ ^(2/3 + /?i)(|m Q | 2 + |m b | 2 ) 2 

+ i(2/3 -/3 1 )(m a -m b ) 2 

+ ^/3 2 (|m Q | 2 Hm b | 2 ) 2 , (8) 

and identify the Ising nematic degree of freedom as 
ip = m a • m b . In the mean-field theory presented in 
this paper, tp is only non-zero in the MS phase. The in- 
clusion of sufficiently strong magnetic fluctuations, how- 
ever, allows the nematic order parameter to be non-zero 
above T/v,— ~— so long as the coefficient of <p 2 in Eq. (|5J) 
is negative, i.e 2/?o < Pi- This is the same condition that 
ensures that the MS phase has lower free energy than the 
OM phase, and hence implies that a nematic transition 
occurs at T > Tn when the SDW state shows stripe or- 
der. It is therefore intriguing that we find that 2/?o — Pi 
changes sign as a function of doping in all cases. This sug- 
gests a strong doping dependence of the nematic phase: 
In particular, for a scenario with a single hole pocket we 
expect that the nematic phase at underdoping will be 
weaker and occur closer to Tjy compared to overdoping, 
or may even be absent. Since this is apparently not ob- 
served experimentally, our results might imply that the 
magnctoclastic coupling plays a more direct role in the 
structural transition^ 



VI. SUMMARY 

In this paper we have presented a weak-coupling study 
of the magnetic order in a two-band model of the iron 
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pnictidcs. Using the Dyson equation for the Green's func- 
tion of an arbitrary commensurate SDW state treated at 
the mean-field level, we have obtained an expansion of the 
free energy valid close to the critical temperature. We 
have shown that this allows three commensurate SDW 
states: The experimentally relevant stripe MS phase, the 
flux OM phase, and the SCO phase for which only one 
sublattice orders. The competition of these phases with 
one another and with IC SDW states has been studied as 
a function of the doping and the variation of key features 
of the non-interacting electronic structure. 

In particular, we have examined systems containing 
two elliptical electron pockets and either a single hole 
pocket at the T point or hole pockets at both the T and 
the M points. In the former case we find that the MS 
state is stable at optimal doping, while a superposition of 
two orthogonal SDW states is stable away from optimal 
doping. In the latter scenario, however, the OM phase 
is realized near optimal doping, while the MS state is 
stable at moderate doping and the SCO state appears at 



higher doping levels. The doping dependence of the asso- 
ciated phase diagrams can be understood in terms of the 
changing nesting properties of the Fermi surface. Our 
results indicate that the single-hole-pocket model is in 
better agreement with experiments, presumably because 
any hole pocket at the M point in the real materials is 
small and poorly nested with the electron pockets. The 
single-hole-pockct picture also suggest that the proposed 
nematic state in the pnictidcs should be highly asymmet- 
ric with respect to electron vs. hole doping. 
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